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Abstract 

We re-analyse the non-standard interaction (NSI) solutions to the solar neutrino 
problem in the light of the latest solar as well as atmospheric neutrino data. The 
latter require oscillations (OSC), while the former do not. Within such a three- 
neutrino framework the solar and atmospheric neutrino sectors are connected not 
only by the neutrino mixing angle ??i 3 constrained by reactor and atmospheric data, 
but also by the flavour-changing (FC) and non-universal (NU) parameters account¬ 
ing for the solar data. Since the NSI solution is energy-independent the spectrum is 
undistorted, so that the global analysis observables are the solar neutrino rates in 
all experiments as well as the Super-Kamiokande day-night measurements. We find 
that the NSI description of solar data is slightly better than that of the OSC solu¬ 
tion and that the allowed NSI regions are determined mainly by the rate analysis. 
By using a few simplified ansdtze for the NSI interactions we explicitly demonstrate 
that the NSI values indicated by the solar data analysis are fully acceptable also for 
the atmospheric data. 

In the appendix we present an updated analysis combining the latest data from 
all solar neutrino experiments with the first results from KamLAND. We show that, 
although NSI still gives an excellent description of the solar data, the inclusion of 
KamLAND excludes at more than Sa the NSI hypothesis as a solution to the solar 
neutrino problem. 
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1 Introduction 


The wealth of data from solar [1,2] and atmospheric neutrinos [3,4,5,6] has 
put neutrino physics in the spotlight and physicists are asking what is it that 
makes solar and atmospheric neutrinos convert. The most popular possibility 
is that of neutrino oscillations, expected in most theories of neutrino mass [7]. 
Indeed an excellent joint description of solar and atmospheric data is obtained 
in this case [8,9]. However, more often than not, models of neutrino masses 
are accompanied by non-standard interactions (NSI) of neutrinos [10]. These 
can have non-universal (NU) or flavour-changing (FC) components, which 
typically co-exist. The simplest NSI does not involve additional interactions 
beyond those mediated by the Standard Model electroweak gauge bosons: it 
is simply the nature of the leptonic charged and neutral current interactions 
which is non-standard because of the complexity of neutrino mixing [7]. On 
the other hand NSI can also be mediated by the exchange of new particles with 
mass at the weak scale such as in some supersymmetric models with R-parity- 
violating [11,12] interactions. Such non-standard flavour-violating physics can 
arise even in the absence of neutrino mass [13,14]. 

Such varieties of non-standard interactions of neutrinos affect their propa¬ 
gation in matter [15]. The magnitude of the effect depends on the inter¬ 
play between conventional mass-induced neutrino oscillation features in mat¬ 
ter and those genuinely engendered by the NSI, which do not require neu¬ 
trino mass [16,17]. These may have a variety of phenomenological implica¬ 
tions [17,18] and have been considered in the context of atmospheric [19,20], 
as well as astrophysical neutrino sources [21]. The impact of non-standard 
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interactions of neutrinos has also been considered from the point of view of 
future experiments involving solar neutrinos [22] as well as the upcoming neu¬ 
trino factories [23,24], 


In this paper we re-consider the case of NSI solutions to the solar neutrino 
problem [25] taking into account the recent charged current measurement at 
the Sudbury Neutrino Observatory (SNO) [2], the 1258-day Super-Kamiokande 
(SK) data and the previous solar neutrino data [1], In contrast to previous work 
we consider a three-neutrino NSI analysis and study also the dependence of 
the solar neutrino NSI solution on the neutrino mixing angle dis whose value 
is presently constrained mainly by the reactor neutrino data [26], Due to its 
energy-independent nature, the NSI solution is in excellent agreement with 
the flat spectral energy distribution observed in the Super-Kamiokande ex¬ 
periment. Thus we focus first on the determination of the allowed solutions 
by considering only the total rates of the solar neutrino experiments. We hnd 
that these solutions provide excellent descriptions of the solar rates, including 
the recent SNO charged current (CC) result. Then we study the impact of 
the day-night data from Super-Kamiokande measurements and show that the 
NSI solutions are consistent with the non-observation of day-night variations. 
Finally, and more important, using simplihed ansdtze for the NSI interactions 
we analyse the impact of the NSI description of solar data on the atmospheric 
data, showing how they are fully acceptable also for the latter. 


This paper is organized as follows. In section 2 we discuss the neutrino evo¬ 
lution and conversion probabilities in the presence of NSI, in section 3 we 
summarize the calculational and £t procedures we adopt. In section 4 we dis¬ 
cuss the impact of the NSI solar solution on the atmospheric data analysis, 
and in section 5 we summarize our results. 
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2 Neutrino Evolution in the presence of NSI 


The most general form of three-neutrino evolution Hamiltonian in the flavor 
base (z/e, t'/i, i^t), in the presence of NSI can be given as 
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where E is the energy, Gp Fermi’s constant, Am ‘^2 = '^3 ~ ''^2 = ^^atm 
7Vg(r) is the number density of electrons along the neutrino trajectory. In our 
calculations of the z/g survival probability we use the electron and neutron 
number densities in the Sun from the SSM [27], while for the Earth matter 
effects we use the density prohle given in the Preliminary Reference Earth 


Model (PREM) [28], 


In Eq. (1) Eij parametrize the deviation from standard neutrino interactions. 
For example is the forward scattering amplitude of the FC 

process z/g-h/ —while \/ 2 GfN represents the difference between 
the non-standard component of the Ue + f and the vi + f elastic forward 
scattering amplitudes. The quantity Nf{f) is the number density of the fermion 
/ along the path r of the neutrinos propagating in the Sun or the Earth. 
We consider two cases, depending on the NSI model, in which the neutrino 
interaction occurs with down-type or up-type quarks. For the case of non¬ 
standard interactions on electrons [22] one would have also to take into account 
the effect of the NSI also in the neutrino detection cross section. 


Note that for simplicity we have set the mass splitting Am^i of the first 
two neutrinos to zero, so that the corresponding mixing angle ^12 can be 
eliminated. This is justihed as we are mainly interested in isolating the effect 
of the NSI in the description of solar neutrino data. In contrast we assume the 
most general NSI flavour structure, consistent with CP conservation. Under 
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these assumptions, we have [29]: 


Cl3 0 Si3 

—-S23'Sl3 C23 S23C13 ) (2) 

~'Sl3C23 ”"523 C23Ci3^ 

where we used the notation Cij = cos'djj and Sij = sin'djj. Similarly to the 
usual oscillation case [30] the Ue survival probability Pee can be written as 

Pee = Cts-Pi* + 4. (3) 

where P®® is the electron survival probability in an effective 2x2 model 
described by the Hamiltonian 



R^VrI = V2 GfN^ 

L J2x2 


V2GFNf{r) 


0 ^eff 


^eS 


with the effective £efr and given in terms of the original Sij parameters as 

= Ci3(ei2C23 — £13523) — 5i3[£23(c23 — ^23) + (^22 ” £33)C23'S23] ; (5) 

= ^22C23 — 2623023^23 + ^33523 + 2Si3Ci3(£i3C23 + £l2523) 

~ ' 5 i 3 (^ 33 C 23 + £22^23 + 2623823023) ■ 

Note that even though we have assumed the most general NSI flavour struc¬ 
ture, the hnal propagation of solar neutrinos can be described effectively as a 
two-dimensional evolution Hamiltonian depending only on the mixing angle 
dis and on two effective NSI parameters (£efr,£efj). 


3 Analysis Method 

In our description of the solar neutrino data [1] we adopt the same analysis 
techniques already described in Refs. [8,31,32] using the latest theoretical stan¬ 
dard solar model (SSM) best-£t fluxes and estimated uncertainties [27]. For 
the neutrino conversion probabilities we use the results calculated numerically 
as indicated above. 
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Since the parameter space is three-dimensional, the allowed regions for a given 
C.L. are dehned as the set of points satisfying the condition 

Xsoi(^efr, 4ff, ^?13) - Xsoi,min < ^X^{C.L., 3 d.o.f.), (7) 

In our numerical calculations we use the survival/conversion probabilities 
of solar electron neutrinos valid over a generous range of NSI parameters 
(eeffj^efr)- the other hand, we use the relevant reaction cross sections and 
efficiencies for the all experiments employed in Ref. [8,31,32]. For the SNO 
case the CC cross section for deuterium was taken from [33]. 

For what concerns the £t of the total rates, we take into account the results of 
the Homestake, GALLEX/GNO, SAGE and Super-Kamiokande experiments, 
together with the latest SNO GG result. Thus we have 5 experimental data, 
which we fit in terms of the 3 parameters Sefr, ^eff and dis; therefore, the total 
number of degrees of freedom is 2. 

In addition to total rates, we also take into account the Super-Kamiokande 
“zenith angle spectrum” [1], which includes both the spectral information 
on the final lepton energy (8 bins) and the information on the zenith angle 
distribution which results from neutrino interactions inside the Earth (1-day 
and 6-night bins). This data sample contains a total of 44 = 6 x 7 -|- 2 bins. 
Thus the number of degrees of freedom is 5 (rates) -|- 44 (zenith-spectrum) - 
1 (free normalization) - 3 (fit parameters) = 45. 

In Figs. 1 and 2 we present the allowed regions in the (Seff, ^eff) parameter 
space for different values of sin^dia, assuming non-standard interactions of 
neutrinos with d-type (Fig. 1) and w-type (Fig. 2) quarks. The shaded areas 
refer to the 90%, 95%, 99% and 99.7% G.L. with 3 degrees of freedom, and 
the best fit point is indicated by a star. Let us now comment some of their 
features. The first thing to notice from Figs. 1 and 2 is that the amount of NU 
is large. Nevertheless this is allowed by experiment and therefore consistent. 
We note also that the NSI strength required in the case of interaction with 
u-type quarks is smaller than for the case of d-type quarks. 

Notice also that can only be large over a very narrow region close to 
0.57, a fact which will be used shortly. More importantly, we note that the 
quality of the solar fit becomes progressively worse as dia increases, in a way 
similar to what happens in a three-neutrino OSG scenario [8,34]. This indicates 
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Figure 1. Allowed NSI regions indicated by the fit of the solar rates (upper 
panels) and global solar data (lower panels) for different values of ??i 3 , assuming 
non-standard interactions of neutrinos with d-type quarks. The shaded areas refer 
to the 90%, 95%, 99% and 99.7% C.L. with 3 degrees of freedom. The best fit point 
is indicated by a star. 
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Figure 2. Same as Fig. 1 but for the case of non-standard interactions of neutrinos 
with u-type quarks. 

that, although with less sensistivity, 'dis can be constrained solely by solar 
data, irrespective of reactor data, also in the context of the NSI mechanism. 
In Fig. 3 we present the dependence of as a function of sin^ 'dia. The solid 
line refers to non-standard interactions of neutrinos with d-type quarks, and 
only the information from total rates is used, while for the dashed line also the 
spectrum-zenith distribution of Super-Kamiokande is included in the £t. The 
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Figure 3. Dependence of as a function of sin^i?i 3 , for different NSI solutions. 


two dotted-dashed lines refer to non-standard interactions of neutrinos with 
M-type quarks, both with and without the spectrum-zenith information. 


Finally in Table 1 we present the best-£t points and goodness-of-ht of oscil¬ 
lation and NSI solutions to the solar neutrino problem. Since the neutrino 
mixing angle dia is strongly constrained by reactor neutrino data [26] we set, 
for definiteness, dia = 0 in what follows. This way we have only two degrees 
of freedom. The numbers in the table refer to a restricted analysis using only 
and the neutrino mixing angle idi 2 for the pure OSC case and only 
and for the pure NSI case. We see that the pure two-parameter NSI 
solution is somewhat better than the corresponding pure OSC solution. 


Solution 

Amii 

tan^(di2) 

^eff 


■ 

Amin 

G.O.F. 

LMA 

2.4 X 10 "® 

0.31 

0 

0 

2.2 

53 % 

SMA 

7.8 X 10 "^^ 

1.6 X 10-3 

0 

0 

4.0 

26 % 

NSI (d) 

0 

0 

3.2 X 10-3 

0.61 

0.60 

90 % 

NSI (u) 

0 

0 

1.3 X 10-3 

0.43 

0.62 

89 % 


Table 1 

Best-fit points and goodness-of-ht of oscillation and NSI solutions to the solar neu¬ 
trino problem, including only rates. 





4 How about atmospheric neutrinos? 


So far we have given a description of solar neutrino data in terms of NSI in¬ 
teractions. We now turn to the issue of atmospheric neutrinos. In principle 
one could imagine a pure NSI description of the atmospheric data itself [19]. 
However, it has recently been shown that, while it may £t the contained at¬ 
mospheric data, such a description can not reconcile them with the up-going 
muon data. The sensitivity of the atmospheric data to the NSI follows mainly 
from the fact that the wide energy range of the atmospheric neutrino data 
sample, from the sub-GeV events up to the up-going samples of MACRO 
and Super-Kamiokande [4,3,5,6] is enough to manifest the energy-dependence 
characteristic of the atmospheric neutrino conversion mechanism. On this ba¬ 
sis it has been shown how, taken altogether, the atmospheric data leave very 
little room for NSI and can tolerate the existence of NSI only at a sub-leading 
level [20]. 


The main question for us then is to evaluate whether the NSI description 
of solar neutrino data is in conflict with the atmospheric data sample. The 
danger lies in the fact that, in contrast to a pure oscillation description of the 
neutrino data [8], even when the value of the third neutrino mixing angle dis 
(presently constrained mainly by the reactor neutrino data) is taken to zero, 
the solar and atmospheric sectors are still coupled with each other through 
the postulated NSI interactions whose strength is fixed, as described above, 
in order to account for the solar data. 


A full description of the atmospheric and solar neutrino data in terms of a 
hybrid OSC-I-NSI description is, at the moment, prohibitive in view of the 
complexity of the problem. It suffices to remind the reader that we expect ten 
relevant independent parameters in this case: in addition to the five relevant 
CP conserving oscillation parameters of Ref. [8] (two mass splittings and three 
angles) there are hve new NSI parameters characterizing the magnitude of FC 
and NU CP conserving non-standard interactions. 


Thus we are forced to make a simplifying ansatz. For definiteness we consider 
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the two following choices of the NSI matrix: 
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With each of these choices the hve new NSI parameters are reduced to the 
two parameters involved in the description of the solar data, as explained 
previously in sections 2 and 3. In fact for both of these choices, once we set 
dis = 0 and '^23 = 45° (see below), we have simply £efr = ^ and = £'. For 
dehniteness we focus here on the case of neutrino non-standard interactions 
with down-type quarks in Eq. (1). 


As for the data here we use the totality of sub- and multi-GeV (e, ju) atmo¬ 
spheric neutrino data [3,4] as well as Super-Kamiokande stop and through- 
going muon data, and MACRO through-going muons [5,6]. 

In the lower panels of Fig. 4 we show the allowed region in the plane (e, Am| 2 ) 
from the analysis of atmospheric neutrino data. The remaining undisplayed 
parameters are set to dia = 0 (indicated by a combination of reactor and 
atmospheric data [26]), i923 = 45° (indicated by the standard atmospheric 
data analysis) and £' = 0.57 (indicated by the solar data analysis presented in 
sections 2 and 3). In the upper panels we show the dependence of as a 

function of e when Am 32 is “integrated out”. 

From Fig. 4 we see that in both cases the value of Am 32 is essentially unaffected 
by the inclusion of non-standard interaction in the > i/g and i/g 

channels, up to values of 10 % or so. This is analogous to what was already 
found in Ref. [20] for the > i/j- channel. We also note that in model (a) 
values of e ~ 0.07 are preferred, since in this case a small contribution of FC 
actually helps in improving the SK contained z/g data. This phenomenon is 
analogous to having a small but hnite solar mass splitting Am^i, see Ref. [35]. 
However, a value of e = 0 is clearly in agreement with the data, as expected 
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Figure 4. Allowed region (90%, 95%, 99%, 99.73% C.L. with 2 d.o.f.) in the plane 
(e, Am| 2 ) from the analysis of atmospheric neutrino data, for each of the two ansdtze 
given in Eqs. (8) (left panels) and (9) (right panels). The best ht point is denoted as 
a star. We fix the undisplayed parameters to 'dis = 0, 1923 = 4:5° and e' = 0.57. The 
top panels show the behavior of as a function of e when Am '^2 is integrated 

out. 

by the excellent quality of the oscillation fit for the pure i Ur vacuum 
oscillations. We see that in both schemes the values of e allowed by the solar 
data analysis are fully acceptable also for the atmospheric data. However in 
model (a) there is a small tension between the solar best £t point (e = 0.003) 
and the atmospheric data (which prefer, with = 1.6, a value e = 0.07), 
so that some slight change in the shape of the solar allowed region might be 
expected when atmospheric data are also included. This effect does not appear 
in ansatz (6), for which any value of £ below 0.03 is practically equivalent. 


5 Conclusions 


In this paper we have re-analysed the non-standard interaction solutions to the 
solar neutrino problem in the light of the recent SNO measurement as well as 
the 1258-day Super-Kamiokande data. In contrast to previous work we have 
used a three-neutrino interpretation of the solar data analysis, displaying how 
results depend on the neutrino mixing angle which is constrained mainly 
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by the reactor neutrino data. More importantly, we have checked consistency 
of the NSI interpretation of solar data with the good agreement that the atmo¬ 
spheric data sample shows with the OSC interpretation. We have found that 
the status of such energy-independent solution to the solar neutrino problem 
is slightly better than that of the OSC solution. While these NSI solutions ex¬ 
ist for reasonable values of the flavour-changing interaction strength ^efr they 
require a somewhat large value of the non-universal parameters, suggesting 
that the solar conversion channel must involve Vr- We have also analysed the 
implications of the solar NSI solution for atmospheric neutrinos, generalizing 
the study in Ref. [20] so as to analyse the sensitivity of atmospheric data also 
to the NSI parameters involved in the solar neutrino channel. By using two 
simplihed ansdtze for the NSI interactions we have explicitly demonstrated 
that the values of e allowed by the solar data analysis are fully acceptable also 
for the atmospheric data. This establishes that the NSI description of the solar 
neutrino data is not in conflict with the atmospheric data sample. For one of 
these models we noted the existence of a slight conflict between the solar best 
fit point (e = 0.003) and the atmospheric data (which prefer, with = 1.6, 
a value e = 0.07), which suggests that some slight change in the shape of 
the solar allowed region might be expected within a fully global description 
of both solar the atmospheric data. Our simplihed ansdtze are justihed to 
the extent that the complexity of the analysis makes a full description close 
to impossible, and possibly un-illuminating. We hud it signihcant, however, 
that the present understanding of solar and atmospheric neutrino data does 
not, as yet, require solar neutrino oscillations or solar neutrino mixing. This 
may have profound implications for model-building, especially in view of the 
difficulties in accomodating bi-large mixing-type solutions within the frame¬ 
work of unihed theories. From an experimental point of view we hud it worth 
pointing that the energy dependence characteristic of the oscillation mecha¬ 
nism has only been demonstrated for the case atmospheric, not solar neutrino 
conversions, as yet. One may argue that our NSI hypothesis is somewhat ar- 
tihcial. However one should bear in mind the fact that most models neutrino 
masses are accompanied by non-standard interactions of neutrinos and there 
are, in fact, some in which NSI are unacompanied by neutrino mass effects. 
The soundness of hybrid schemes such as the ones indicated here may indi¬ 
cate more subtle new ways for accounting for the presently observed neutrino 
anomalies from hrst principles. We look forward to new solar neutrino exper¬ 
iments to probe neutrino NSI with improved sensitivities, such as suggested 
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in Ref. [22], as well as to the prospects of future neutrino factories performing 
correlated OSC-NSI studies, as suggested in [23,24], 
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Appendix: Implications of the KamLAND result. 


The original version of the present paper showed that a pure NSI description 
of solar data with massless and unmixed neutrinos is slightly better than the 
favoured LMA-MSW solution. Moreover, the NSI values indicated by the solar 
data analysis do not spoil the succesful oscillation description of the atmo¬ 
spheric data, establishing the overall consistency of a hybrid scheme in which 
only atmospheric data are explained in terms of neutrino oscillations. However, 
the recent results of the KamLAND collaboration [36] reject non-oscillation 
solutions in the solar sector under the assumption of CPT invariance, in such 
a way that solutions based on NSI should be strongly disfavoured by Kam¬ 
LAND results. In this appendix, we re-evaluate the status of NSI solutions of 
the solar neutrino problem, in the light of the first KamLAND data. 

The KamLAND experiment is a reactor neutrino experiment whose detector 
is located at the Kamiokande site. Most of the W flux incident at KamLAND 
comes from nuclear plants at distances 80-350 km from the detector, making 
the average baseline of about 180 km, long enough to provide a sensitive probe 
of the LMA-MSW region. The target for the W flux consists of a spherical 
transparent balloon hlled with 1000 tons of non-doped liquid scintillator, and 
the antineutrinos are detected via the inverse neutron /3-decay process W + 
p ^ + n. The KamLAND collaboration has for the first time measured 

the disappearance of neutrinos produced in a power reactor. They observe 
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a strong evidence for the disappearance of neutrinos during their flight over 
such distances, giving the first terrestrial conhrmation of the solar neutrino 
anomaly and also establishing the oscillation hypothesis with man-produced 
neutrinos. 

In this appendix we perform a combined analysis of the solar neutrino data 
with the hrst KamLAND results in terms of neutrino non-standard inter¬ 
actions with matter, setting for dehniteness dis = 0. The analysis of solar 
neutrino data includes the most recent results from all experiments [37], as 
well as the latest measurements from SNO [38] presented in the form of 34 
data bins (See Ref. [39] for a detailed explanation of the new solar analysis). 
The details of the theoretical Monte-Carlo and statistical analysis of the Kam¬ 
LAND results are given in Ref. [40]; in particular, the KamLAND y^-function 
is calculated assuming a Poisson distribution of the experimental data, as 
described in Sec. IV of that paper. 

In Table 2 we present the best-£t points for the oscillation and NSI solutions 
to the solar neutrino problem obtained in our analysis, together with its cor¬ 
responding value of The value of confirms our previous result: NSI 
picture gives a good description of the solar neutrino data. In fact, for the 
case of neutrino NSI with up quarks we obtain a slightly better fit than for 
the LMA-MSW solution. 

The inclusion of the latest solar data sample [37,38] does not change the status 
nor the allowed regions of the NSI solution in any significant way. In contrast, 
the inclusion of KamLAND data in the analysis changes dramatically the situ¬ 
ation. Note that for the KamLAND experiment matter effects are very small, 
and can therefore safely be neglected. As a result, the NSI solutions predict 
no reduction in the z/g flux, in conflict with what is observed at KamLAND. 
This worsens the status of the NSI hypothesis, (reflected in the corresponding 
Xsoi+KL values, in Table 2) with respect to that of the LMA-MSW solution, 
which predicts the correct suppression factor observed at KamLAND. From 
our analysis we obtain that the description in terms of neutrino NSI with down 
quarks is rejected with a Ax^ = 14.2, corresponding to 3.8(T, with respect to 
LMA-MSW, while NSI of neutrinos with quarks of type up are rejected at 3.6a 
level (Ax^ = 12.8) as the explanation for the solar neutrino anomaly plus the 
KamLAND disappearance of neutrinos. 
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Solution 

l> 

tan ^ ( 1912 ) 

^eff 

^eff 

Xsol 

aLi+kl 

LMA-MSW 

7.2 X 10 "^ 

0.46 

0 

0 

65.8 

71.9 

NSI (d) 

0 

0 

3.2 X 10"3 

0.62 

66.9 

86.1 

NSI (u) 

0 

0 

1.3 X 10-3 

0.44 

65.5 

84.7 


Table 2 

Best-fit points of oscillation and NSI solutions to the solar neutrino problem before 
and after the inclusion of the KamLAND data in the analysis 
In summary, we have shown that the inclusion of KamLAND data in the 
analysis of the solar neutrino anomaly excludes the interpretation based on 
neutrino non-standard interactions with matter at more than 3a. Therefore, 
non-standard interactions may at best play a sub-leading role in solar neu¬ 
trino propagation. In fact, one may use the conhrmation of the LMA-MSW 
oscillation solution together with the experimental data from KamLAND and 
solar neutrino experiments in order to determine improved restrictions on NSI 
parameters. 
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